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negative electrode. Using nitrogen data for a numeric example, the following observations 
are made: a sufficiently strong applied electric field pushes the ion density toward that of the 
electrons through a well - a constrictive phenomenon. Both a transition region, dominated 
by density gradients, and a diffusion-driven zone are found to move the system toward the 
plasma interface. The characteristics of this process are influenced by the applied electric 
field, but the instability of the chemistry-induced stagnation regions precludes numeric 
convergence. Insufficient dissipation may prevent the stability of the anode fall model as 
presented. Suggested improvements to the model descriptions include considering the effects 
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L INTRODUCTION 



Whenever an object starts moving in a given direction, a second object must be 
present which accepts the reaction of the force that accelerates the first, so says Newton’s 
third law of motion. When a spacecraft moves through empty space, the accelerating forces 
act between the vehicle and the mass of the propellant expelled. In 1895, the distinguished 
Russian rocket pioneer Konstantin Eduarovitch Tsiolkovskii derived the famous “Tsiolkovskii 
equation,” the basis of all theoretical work on rocket propulsion. Relating the burnout 
velocity of a rocket with its exhaust velocity and its mass ratio, this “rocket equation” shows 
how important it is for the exhaust velocity of a rocket motor to be as high as possible. [Ref. 
1. pp. 2-3] 

Since mankind’s advent into space, several types of space flight propulsion systems 
have been developed to include chemical, nuclear, electric and solar propulsion devices with 
the majority of spacecraft to date employing chemical rockets. While the science and 
technology of chemical rockets developed at a rapid pace, the idea of using the benefits of 
Maxwell’s relations to eject charged particles from a thrust chamber progressed slowly with 
long intervals between periods of creative thought. In 1945, Herbert Radd recommended the 
use of fast electrically accelerated particles to reduce fuel mass and in 1946, Jakob Ackeret 
analyzed the performance capability of a spacecraft by arguing that the total kinetic energy 
carried away by the exhaust beam corresponded to a reduction of the mass of the vehicle 
according to E=mc 2 . Hence, assuming that a certain fraction of the vehicle mass is converted 
into kinetic beam energy, the total propellant mass can be calculated which in turn leads to 
a maximum terminal velocity of the spacecraft, a result which is independent of particle size, 
time of propulsion, conversion method and magnitude of the vehicle. In essence, the 
optimum exhaust velocity is determined simply by maximizing the terminal velocity. [Ref. 1 : 
pp. 2-3] 

Propulsive devices that use Maxwell’s Law’s differ from chemical systems: the latter 
is characterized by having the same energy and propulsive sources, while in the former the 
energy conversion method differs from that of propulsion. In essence, electric propulsion 
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systems allow greater carrying capacity with shorter flight times. For example, the Hohmann 
transfer time from a Low-Earth-Orbit (LEO) to the orbit of Pluto requires approximately 45 
years if chemical means are considered; an electrically propelled vehicle, in comparison, 
requires approximately three years. [Ref. 1: p. viii]. The large difference in mass payload ratio 
(final mass/initial launch mass) obtained is highlighted by an analysis for a Mars mission: a 
chemical system using a high impulse Hohmann trajectory from LEO delivers approximately 
10% to 18% of the launch mass to the Martian surface [Ref. 2: pp. 152], while an electric 
system using a low impulse spiral trajectory could deliver 20% to 60% of the launch mass, 
depending on the desired transit time. With effective exhaust velocities as high as 10,000 m/s, 
electric propulsion thus offers the performance envelope needed for manned interplanetary 
missions. [Ref. 2: pp. 37] 

Electric propulsion systems are divided into three categories: 1) electrostatic, where 
acceleration is achieved through the interaction of coulombic fields with charged propellant 
particles, 2) electrothermal, where propellant heated electrically is expanded 
thermodynamically through a nozzle, and 3) electromagnetic, where acceleration is achieved 
through the interaction of the Lorentz force with highly ionized gases (plasma). Of the types 
of electric propulsive thrusters available, electrothermal devices include the resistojet and 
arcjet. It should be noted that the current necessary to form the plasma in these devices also 
allows the magnetic field to accelerate the gas. These systems (and for that matter, any 
device using a thermal nozzle) are, however, limited by the amount of heat that can be added 
to the flow and are thus ultimately dependent on the materials used. As electromagnetic 
(EM) thrusters rely on Lorentz interactions to accelerate the flow, these devices are not 
limited in the sense of electrothermal systems and therefore offer different operational 
characteristics. The various types of EM thrusters available include the ion engine and plasma 
thrusters (wide variety to include electrodynamic, magnetoplasma-dynamic or MPD, pulsed- 
plasma. Hall accelerator, Lorentz force accelerator, pulsed coaxial thruster, plasma arcjet, 
pulsed ablative thrusters, to name just a few). [Ref. 2: pp. 565-596] The MPD thruster is 
relevant to this work. 

MPD systems are a relatively recent addition to the electric propulsion family. 
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originally discovered when peculiar behavior of an electrothermal arcjet thruster was 
investigated. Researchers found that without significant mass flow-rate applied, the arcjet 
continued to operate in a mode clearly unique: pressure within the thrust chamber was too 
low for the device to be choked, and so it was theorized that an induced magnetic field from 
the high arc current accelerated the plasma flow - the magnetoplasmadynamic thruster was 
“bom.” [Ref. 3] 

The area of research to which MPD thrusters are applied is termed 
magnetohydrodynamics (MHD), the field of physics that investigates how moving conductive 
(electrolytic) fluids behave while subjected to applied electric or magnetic fields. The fluids 
involved must be conductive, or ionized sufficiently to pass a current, and in the case of an 
MPD, the fluid must reach a plasma state. MHD devices behave analogously to motors or 
generators and can either add or remove energy from the fluid in exchange for electrical 
energy; systems developed include the MHD tunnel drive for submarines and the MHD 
turbine. Magnetoplasmadynamic (MPD) thrusters are currently being considered for 
application to primary propulsion system on robotic and piloted interplanetary missions; the 
high specific impulses offered by MPD devices can reduce launch mass requirements by over 
a factor of two compared to chemical system if the power to thrust conversion efficiency is 
over 50%. [Ref. 4] 

MPD thrusters are classified by considering two major features: 1) the power supply 
which defines the continuous, quasi-steady and pulsed MPD systems, and 2) the magnetic 
field which characterizes the self-field and applied-field MPD thrusters [Ref. 3], A 
continuous, or steady-state, thruster is one in which the propellant flow and current 
throughout the device are continuous, thereby allowing the system to reach a high 
temperature equilibrium state. The quasi-steady thruster does not reach equilibria since it 
uses a continuous gas flow with rapidly pulsed current, through which high amounts of energy 
are provided in short bursts. The pulsed thruster is similar to the quasi-steady MPD system 
except it uses less frequent pulses; equilibrium flow is also not reached and consequently low 
temperature materials may be used for constructing such a device. The applied- and self-field 
thrusters define the way in which the primary accelerating field (magnetic) is affected. The 
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applied-field system employs a large current in an external coil, or large permanent magnets, 
to accelerate the plasma whereas in the self-field device the current not only ionizes the gas 
but also creates the field that accelerates the resulting plasma. 

The magnetoplasmadynamic thruster was originally discovered in the 1960's, and yet 
much work remains to understand its fundamental operation. For example, the steady state, 
self-field MPD thruster operated by the University of Southampton, U.K., was originally 
developed to simulate the atmosphere effects on structures in LEO, an environment 
dominated by highly reactive atomic oxygen. The MPD is suited to simulate such conditions 
as it offers a unique combination of high thrust and high exit velocity. The Southampton 
MPD presently operates at power levels in the range 6-17 kW but suffers from electrode 
erosion, instabilities, a contaminated beam and short run-times. Work on this thruster has 
identified electrode erosion as a major problem, making the device unsuitable for its intended 
function. Specifically, the mode of erosion differs from that seen at other institutions and the 
levels of erosion are higher and exhibit a rate increase not seen with other devices. [Ref. 5] 
Electrode erosion, current spotting, frozen flow loses and electrode power deposition are just 
several factors that limit the performance of MHD devices [Ref. 6, 7]. Specifically, anode 
power deposition is the single largest power loss mechanism in MPD thrusters operating at 
sub-megawatt power levels. [Ref. 8, 9, 10, 11, 12] 

To gain further insight to MPD phenomena, computer codes that accurately describe 
observed data from steady-state MPD thrusters have been developed [Ref. 13, 14, 15], 
However, these codes do not adequately describe observed data from quasi-steady thruster 
experiments. It has been suggested that the lack of proper electrode modeling (i.e., sheaths 
and fall potentials) in these codes may explain this discrepancy [Ref. 16]; perhaps the difficult 
set of coupled, nonlinear partial differential equations involved has limited analytical work to 
model the anode sheath and ambipolar regions. For example, Hugel [Ref. 17] and 
Subramanian [Ref. 18] address the influence of the sheath region, but do not model the 
electric field, temperature or sheath fall voltage. Hence, the focus of this work: analysis of 
the anode voltage drop region. 
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II. LITERATURE REVIEW 



The principal voltage loss mechanism near an electrode can be divided into two main 
classes: ohmic and sheath losses. Ohmic resistivity occurs because of the finite conductivity 
of areal plasma, particularly thermal boundary layers, degree and kinetics of ionization and 
Joule heating. Sheath drops are caused by Debye shielding which forms a space charge field 
adjacent to the electrode [Ref 19: pp. 200-204]. Furthermore, material problems restrict the 
temperatures at which the electrodes can operate. In many cases cooling of the electrodes 
is required since the plasma, in order to maintain a high ionization, must be hotter than the 
working temperatures of most materials. This temperature difference between the electrodes 
and the plasma further aggravates the voltage losses because of the presence of the thermal 
boundary layers. As has been shown, voltage losses can exceed 50% of the total power 
output, with sheath drops accounting for a large fraction of this loss. [Ref. 6, 7] 

The tool most commonly used for measurement of local plasma properties is the 
electrostatic probe, a small electrode which may be biased positively or negatively with 
respect to the plasma in which it is immersed. Such work is relevant to non-emitting anodes 
which, though heavily biased, also probe and disturb the plasma locally. As a result, most 
analytical work on sheath phenomena is embodied in probe theory investigations where the 
effects of an electrode (probe) locally disturb a quiescent plasma. [Ref. 20,21: pp. 1-7, pp. 
8-9] 

In 1924 Langmuir pioneered the use of electric probes, giving the first sound 
mathematical basis to this diagnostic technique. The analysis of electrostatic probes in 
collision-dominated (higher pressure) plasmas began with Davydov and Zmanovskaja [Ref. 
22] who assumed that a quasi-neutral diffusion controlled region extended to within one mean 
free path of the probe where it was matched to the edge of a free-fall sheath; a transition 
region between the two regimes was not attempted, however. Ecker, et al. [Ref. 23] 
attempted matching the solution in the collision-dominated region directly to that in the free- 
fall sheath using variational principles, work which was criticized by Cohen [Ref. 24] based 
on mathematical grounds: something was forced to fit via some artificially imposed boundaiy; 
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matching between separate regions was accomplished through variational principles and 
consequently was not asymptotic in nature [Ref. 24], 

Existing sheath solutions are limited to special case simplifications. Several solutions 
are available for spherical probes by Kiel [Ref. 25], whose analysis neglected diffusion, Barad 
and Cohen [Ref 26], Su and Lam [Ref. 27], as well as by Cohen [Ref. 24, 28], The latter 
group carried out the first completely systematic analyses of spherical probes in collision- 
dominated plasmas. These researchers assumed the continuum-equations valid everywhere 
in the plasma right up to the probe surface; their results showed both the sheath and the quasi- 
neutral regions appeared naturally from the diffusion equations. In these efforts, two major 
assumptions prevail: the plasma was assumed so lightly ionized that only charge-neutral 
interactions were accounted, and the continuum equations were postulated with an isothermal 
formulation. For cool, nonequilibrium plasmas, however, charge-charge collisions must be 
considered [Ref. 29: p. 165]. To that end, McKee and Mitchner [Ref. 30] investigated a 
collisionless sheath with constant transport properties that also included ionization as well as 
recombination effects in the ambipolar region. Exploring spherical probe theory further, 
Cohen and Schweitzer [Ref. 31] then extended Cohen’s original asymptotic analysis [Ref. 24] 
for the entire region, an analysis which allowed for weak ionization and recombination effects 
while also assuming constant transport coefficients. Su and Sonin [Ref. 32] next formulated 
a spherical probe theory by considering the effects of electron-ion collisions by approximating 
diffusion coefficients as a function of the neutral particle density, a density assumed as 
constant. Barad and Cohen [Ref. 26] then presented a continuum theory (for spherical 
probes) in which non-constant transport coefficients were used, coefficients based on charge- 
charge collisional assumptions. This analysis, however, did not incorporate ionization or 
recombination effects. In work similar to Lam’s original theory for flow of weakly ionized 
gases [Ref 28], Stahl and Su [Ref. 33] use the same approach of separating the sheath, 
ambipolar region and free stream; their methods prove the existence of a sheath on a flat 
probe, yet plasma chemistry is not considered (ionization and recombination). It should be 
noted that both Su and Lam [Ref. 27], for spherical probes, and Stahl and Su [Ref. 33], for 
flat probes, consider transition regions to match the inner sheath to the outer quasi-neutral 
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or ambipolar zones - albeit both research groups neglected plasma reactivity in their 
formulations. In related work, Godyak and Sternberg [Ref. 34] applied the ionization 
potential as the boundary condition for both the plasma and the sheath in an effort to 
understand the electrodynamic properties of the plasma-sheath boundary for the cathode; their 
analysis, however, assumed constant ionization for the plasma reactivity with recombination 
neglected. 

It may be inferred from these efforts that anode-voltage-drop-region analyses remain 
insoluble due not only to the nonlinear nature of the equations, but also to the transitional 
nature of particle motion in most sheaths. In an extensive review of probe theory, Chung, 
Talbot and Touryan [Ref. 20] state that a general solution for determining charge density and 
species temperature for probes small relative to boundary layer thicknesses is not available. 
Additionally, Nasser [Ref. 35] discusses an elementary theoretical approach to the glow 
discharge problem by introducing a set of one-dimensional governing equations (continuum) - 
the electron and ion current equations, the net production of electrons and ions, and the 
Poisson equation. According to Nasser, these equations should have a solution, but most 
attempts to achieve such have failed, with the difficulty lying in the boundary conditions [Ref. 
35: p. 412]. Voltage losses at electrode boundaries and surface erosion, together with sheath 
effects which influence ionization, play an important role in plasma devices and must be 
controlled in designs of practical interest. Particularly, thermal arcjets and MPD accelerators 
deposit between 15% and 80% of the input power into the anode. This presents not only a 
severe performance penalty, but also introduces thermal design problems since the heat thus 
generated must be radiated away from the thruster surfaces [Ref. 7, 9], 

Completely satisfactory solutions to the anode voltage drop phenomena do not exist 
for several plasma conditions; one such case involves steady, collisional, low-temperature 
isothermal discharges [Ref. 36]. For instance, the effect of several parameters on the region 
require further study, including the influence of temperature, the role of the problem 
dimensionality as well as the effect of plasma chemistry. It remains unclear, moreover, as to 
the exact cause of the anode spots which have been observed for various MPD thruster 
operating conditions, a behavior which appears roughly as periodic oscillations [Ref. 6, 7, 9], 
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A one-dimensional, collisional, equilibrium solution can satisfactorily reproduce the observed 
electric field and charge density distributions for the entire anode voltage drop region where 
the electron temperature equals that of the heavy species [Ref 36], However, this model 
cannot describe any decrease in current density away from the surface, or current constriction, 
at the anode surface which might be necessary in non-equilibrium. A two-dimensional model, 
developed by Biblarz and Dolson [Ref 6], represents these phenomena, in the absence of 
species ionization and recombination, and predicts the voltage drop in the region: the sheath 
accounted for a majority of the anode voltage drop [Refs. 6, 7, 36], 

Since the pioneering work of Tonks and Langmuir (relating to low-pressure 
discharges), the anode description has remained as one of the classic problems in plasma 
physics. In this work, the high density problem is generalized to any electrode surface in 
contact with a partially-ionized plasma when the current flow is sufficiently small: the 
problem of a non-emitting electrode (probe) in contact with a thermally generated plasma. 
An anode is explicitly considered along with a non-emitting cathode and the governing one- 
dimensional continuum equations are formulated. The electrodynamic behavior of plasma 
reactivity in the electrode regions is then discussed with considerations given to three distinct 
zones of the anode voltage drop regime: inner or space-charged sheath, transition and outer 
or quasi-neutral ambipolar zones. 
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III. GOVERNING EQUATIONS AND NUMERICAL DATA 



A. GOVERNING EQUATIONS 

A flat electrostatic probe immersed in a high-density, flowing, partially-ionized gas is 
considered. For plasmas of moderately high density and low degree of ionization (as applied 
to discharge lasers and MUD devices), the plasma region affected by an electrode, i.e., the 
electrode region, can be considered of the boundary-layer type with the transport of charges 
approximated as one-dimensional. Traditionally, the electrode region is further separated into 
two distinct zones, as delineated by the species densities: a sheath and an ambipolar region, 
both which may be of the order of tenths or even hundredths of a millimeter in extent [Ref. 
6] where the sheath extent is driven by variation of the electric field and the quasi-neutral 
region extent is influenced by the influences of chemistry. 

The problem is formulated by considering the hydrodynamic continuum equations for 
two species, ions and electrons. Induced and applied magnetic fields are neglected with 
energy conservation implicitly satisfied since heat transfer, radiation. Joule heating and surface 
emission effects are at first neglected (for example, B=0, VT e =0, VTj=0, J-E=0). 

The equations below are assumed valid throughout the entire domain of interest, to 
include the electrode (probe) surface. Table (1) lists the nomenclature used: 



d ( Mn . v . ) 

dt 

a ( mn e v e ) 

dt 



+Mn . v.'Vv. = qn . E-kT.Vn . -Mn .f.v.-F. 
111^-1 11 111 1 


Equation (la) 


+ mn v Vv = -qn E-kTVn -mn f v -F 

eee^e ee eeee 


Equation (lb) 


V\E= _SL ( n . -n ) 
e 1 e 

O 


Equation (2) 


Vj. = qV(n.v.) 
Vj e = qV(n e v e ) 


Equation (3 a) 
Equation (3 b) 


V ( n. v . ) =vn e -a 2 n.n e 
V(n v )=-vn +a,n,n 


Equation (4a) 
Equation (4b) 


V* J = 0 J=j+j 


Equation (5) 
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Equations (la) and (lb) describe the momentum-transfer between the ions and electrons with 
neutrals: the left-hand side represents the unsteady and convective contributions, 

respectively, while the right-hand side symbolizes coulombic, pressure, collisional and 
frictional forces, respectively. Next, the electrode is considered as a positive, non-emitting 
surface (anode). Through this process, a flow of negative charges is induced toward the plate 
and a flow of positive charges is repelled away from the plate; the anode, in completing the 
arc current with the cathode, collects incoming electrons while driving all positive charges 
toward the negative surface. This separation of charge in the anode fall region produces a 
divergence of the electric field and a variation of potential away from the plate, as given by 
equation (2) (the variation of potential is readily deduced from Gauss’s equation (2)). The 
current densities are denoted by equations (3a) and (3b) in terms of the species flux of 
charged particles per unit area. Equations (4a) and (4b) are then the continuity equations for 
the ion and electron flux as caused by the ionization and electron-ion (two-body) 
recombination processes: positive ions are produced within the anode fall region by electron 
impact with neutrals and conversely, positive ions combine with negative electrons to form 
neutrals. Equation (5) is the conservation of total charge; the sum of the species currents 
remains constant for a one-dimensional formulation. Hence, from set (1-5), the anode model 
presented in this work is developed with the following assumptions: 

1 . Steady state; free-stream plasma is neutral with uniform charge densities (Saha 
equilibrium). 

2. Mean free path of the charged particles is small compared to the sheath thickness 
which in turn is much smaller than the dimensions of the body; the sheath is treated 
as a continuum. 

3. Chemical equilibrium in the boundary layer; the two-body recombination 
coefficient is held constant over the entire anode fall region as a first approximation 
[Ref 29: p. 238]. Further, the ionization coefficient is treated as Townsend’s first 
coefficient. 

4. There is no applied magnetic field and induced magnetic field effects are neglected 
(low Magnetic Reynolds number) - consequently no ion slip. 
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5. No continuum radiation losses, no ion emission, no significant Joule heating 



6. Isothermal: VT e =VT£=VT o =0 and T e =T i =T 0 . The true electron temperature will 
either be close to the anode surface temperature or the free stream temperature since 
thermal equilibrium is assumed in the boundary layer. Because of the relatively small 
number of collisions that the electrons suffer in the boundary layer, it is frequently 
argued that the electron temperature is close to the free stream temperature [Ref. 37], 

7. Difiiision velocities of the ions and electrons due to electric field interactions are 
small in comparison with their thermal velocities; diffusion coefficients are assumed 
constant. 

8. End effects, in the case of a planar surface, are neglected. 

9. Under the assumption of a highly biased probe, it is possible to ignore convection 
not only in the very thin transition region, but also in the sheath region, even for 
sheaths with the same order of thickness as the boundary layer [Ref. 33], Assuming 
perfectly elastic collision between the species and neutrals, essentially a constant 
collisional frequency is assumed sufficiently large so that the convective derivative in 
(1) can be neglected. 

10. With the assumptions of steady state and negligible convection, a spatially 
homogenous fluid results. 

For a semi-infinite plane with coordinate outward from the planar positive surface so 
that V=d/dy, set (1-5) thus becomes: 
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Here, equations (3a) and (3b) have been incorporated with equations (4a) and (4b). Using 



species mobilities, a = — 2 — and the diffusion coefficients, D.= — — (the Einstein 

’ m iAe ' m iAe 




flux, r ie , where Fick’s law of diffusion is a special case (E=0 so that g^ c =0). Specifically, (1) 
and (2) become 



regions occurs simply because more particles begin their motion in the more dense area - the 
flux is proportional to the density gradient. [Ref. 38: pp. 135-138.] But as the species travel 
toward the electrode, reactivity occurs - ions and electrons recombine while electron collision 
with neutrals breaks neutrality adding charges to the electrode region; flow neutrality is 
broken and the fluid is no longer a plasma. Thus, under constant-property conditions, the 
species continuity equations together with Gauss’s equation describe the flux of charges in 
the vicinity of an electrode or plasma probe. [Ref. 33, 39] With the application of the Einstein 
relation, set (1-5) are re-cast to give the model used in this work, the geometry for which is 
shown in Figure (1): 




Essentially, diffusion is a random process in which a net flux from dense regions to less dense 
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Figure 1 : Geometry of Anode Fall Region 
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To obtain the mathematical model used in this work, the definition of the diffusion coefficient 
is used to re-write the species elastic collision frequencies, ^[kTJ/pDjinJ for isothermal 
conditions, and the current density is defined as j s =qi\v s where s=i, e. The resulting equations 
(1-5) are then solved for their highest derivative terms: 
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Equation (7b) 
Equation (8) 

Equation (9a) 
Equation (9b) 



In the formulation of equations (7-9), 14 variables result (three universal constants 
included). Before analysis of this system, the equations must be written in terms of 
dimensionless variables and proper parameters. For this work, both the Buckingham Pi 
Theorem and Fractional Analysis techniques were applied (Appendix A). The Buckingham 
method, although a useful tool for revealing the dimensionality of the problem, does not bring 
forth the physical effects which influence system behavior. On the other hand, fractional 
analysis coupled with the use of the governing equations, begins to give physical meaning to 
the dimensionless parameters obtained. 

Using undisturbed plasma values for normalizing the variables gives three primary 



scales for the 
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the length scales listed in Table 3: 



Governing Equations 
[Equations (7-9)] 
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Table 2: Non-dimensional Parameters 
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Table 3: Length Scales via Fractional Analysis 
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It should be noted that the length scales 1 ^, L n and Lj are in agreement with reference 50 if 
the species conduction currents in the undisturbed plasma are written in terms of thermal 
velocities and the above electric field and number density scales are combined to form the 
Debye shielding length for planar surfaces ( L E L n =A D 1/2 ). 

B. NUMERICAL DATA 

To gain further insight into the physical models presented, numerical analyses were 
conducted with data from reference 36 for 6000°K nitrogen; data for the ionization and two- 
body recombination coefficients are found in Appendix B. 





Data Value 


Temperature (°K), T 0 


6.0000e+03 


Total Current (mA/cm 2 ), J: used for j^, 


1.0000e+02 


Undisturbed Plasma Charge Density (1/m 3 ), o^: used for n^, n^ 


1.0000e+19 


Undisturbed Plasma Total Density (1/m 3 ), N: 


2.0000e+24 


Undisturbed Plasma Electric Field (V/m), E„ 


1.2000eK)4 


E/N (Vcm 2 ) 


6.0000e-17 


Pressure (N/m 5 ), p: from [p=NkT] 


1.6560e+05 


Two-body recombination coefficient (m 3 /s), a 2 : [Ref. 40] 


1.0000e-13 


Ionization coefficient (1/s), v: from [Ref. 41 ] with v/N= 10' 18 (m 3 /s) for above 
E/N converted to Td (lTd=10‘ J1 Vm J ) 


1.0000e+06 


Ion Diffusion Coefficient (mVs): D,N=4.46(10 18 )T o 1/2 from [Ref. 42: pp. 2-28] 


1.7274e-04 


Electron Diffusion Coefficient (m 2 /s): D/D e ~ni/M 


3.1665e-01 



Table 4: Numerical Data for Analysis - Nitrogen 
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IV. PHYSICAL MODEL 



A. SPACE-CHARGE REGION: SHEATH 

Voltage losses at electrode boundaries, surface erosion and sheath effects play an 
important role in plasma devices and must be controlled in designs of practical interest. 
Particularly, thermal arcjets and MPD accelerators deposit between 15% and 80% of the 
input power into the anode. This presents not only a severe performance penalty, but also 
introduces thermal design problems since the heat thus generated must be radiated away from 
the thruster surfaces. [Ref. 7, 9] So a natural question: what happens at the anode? 

Near a surface (anode) in contact with a plasma there is a region, the sheath, in which 
the electric field is strong and in which the electron number density is much larger than that 
of the ions (see Figures (2) and (3)). 




Figure 2: Axial Potential Profile for a High 
Current Arc Between Two Electrodes (not 
to scale) [after Ref. 3: p. 1 14] 




Distribution of Species Densities - 
Anode Fall Region (not to scale) 



Figure (3) shows the qualitative behavior, in the anode fall region, of the spatial 
distribution of species densities; it is expected that from the wall outward, the densities 
approach one-another, eventually reaching Saha equilibrium at the plasma boundary. The 
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overall shape of the curves (Figure (3)) stem purely from the density flux — = \n- a 2 n 2 , 

dt 

in the absence of any other types of influences. However, it is usually the interaction of a 
systems’ many parts, not the action of each piece independently, that determines the final 
system behavior. 

The mathematical model developed for the sheath, whether anode or cathode, depends 
on the details desired: whether surface effects are considered, whether induced or external 
magnetic fields are included, as to how the flow chemistry is modeled, among others. For the 
cathode, sheath models usually adopt the assumptions that electrons are in equilibrium with 
the electric field and that the ionization in the sheath can be neglected [Ref. 34], In 
comparison, the electrons in the anode fall must supply by ionization enough ions to account 
for the ion current that flows out of the positive column toward the cathode [Ref. 43: pp. 59- 
60]; electrons emerge from the column, are attracted by the anode, gain energy through the 
voltage drop and ionize neutrals. The extent of the anode sheath is then determined by the 
space charge distribution and by the magnitude of the gas ionization potential [Ref. 35: p. 
420]; hence, the associated length scale for the sheath is governed by that which results from 
the dimensional analysis of the Gauss equation (8). 

Moreover, after a current flow is established between the plates, the anode region may 
exist in a vapor that issues from the electrodes. In vacuum arcs. Miller [Ref. 44] characterizes 
the anode region as operating in one of five distinct modes, ranging from a passive, low 
current mode to a high current, fully developed spot mode. These observed instability 
phenomena are similar to those noted in MPD thrusters - anode spot formation at high current 
clearly limits anode lifetime [Ref. 6, 7, 9, 17]. What could cause such “spotting?” 

In a comprehensive review of anode-spot phenomena. Miller suggests [Ref. 44] that 
a sudden increase in the voltage noise frequently indicates that an anode foot-print could be 
forming or that a transition into the anode-spot mode could be occurring. Further, this abrupt 
change in voltage has been associated with a sudden change of ion density in the region near 
the anode, a density decrease that would leave the local electron space charge uncompensated 
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and thus produce a low-conductivity region. The sudden increase in voltage would then 
reflect this local shortage of ions - an ion depletion region or ion starvation region. Even- 
though the most suitable model for the anode region is still actively debated, the general idea 
of the appearance of an ion deficiency region near the anode as triggering the transition into 
the anode-spot mode has been fairly supported. [Ref. 44] So, to better understand the 
processes involved within the anode fall, equations (7-9) are modified as follows: 

Because the anode-spot instability is likely related to the species density distribution 
within the anode fall region, the role of net charge production, the chemistry, should be 
critical to the overall system behavior: from equations (7-9), the species densities are 
influenced by both coulombic interactions as well as by the currents and consequent 
production phenomena, the fine balance between ionization and recombination. Therefore, 
the current densities as given by equations (9a-9b) are combined with the species densities, 
equations (7a-7b): the latter are differentiated with respect to the spatial variable in this effort. 
The result gives for the space-charge region a system of two second order, ordinary, non- 
linear, coupled differential equations with one first order, coupled, ordinary differential 
equation. The non-dimensional parameters given in the bracket notation are listed in Tables 
2 and 10: 
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Essentially, the conservation laws are combined: species momentum now reflects the 
influence of charge production terms. Here, the left-hand side of equations (10a- 10b) is the 
rate of change of the number density gradient, an accelerating behavior. The gradient of the 
species number density is found in the first term on the right-hand side of set (10a- 10b), a rate 
influenced by the strength of the electric field as measured by the imbalance between the 
species number densities, the Gauss equation (10c). The first term of the right-hand side of 
(10a- 10b) describes the interaction of coulombic and gradient forces, the influence of the 
equation of motion as described in the earlier formulation of (1-5); the third term completes 
the formulation of motion in this model. The general concept of energy conservation begins 
to take shape with the second and fourth terms of (10a- 10b): The former has included two- 
body recombination effects as diminished not only by coulombic interactions, but also by 
diffusion. The fourth term illustrates a source of energy to the system , ionization, again as 
influenced by diffusion. Essentially, the accelerating terms on the left-hand side of (10a) and 
(10b) are described by balances of motion (gradients) and energy (ionization, recombination) 
where the exchange of energies is brought about and influenced by diffusion processes. It 
should stand to reason that as the species density gradient increases, energy losses should 
become more dominant until the gradients diminish their contribution to the flow: 
recombination should begin to take hold until eventually equilibrium is restored, balancing 
diffusion with energy source and sink. 

B. TRANSITION ZONE: QUASI-NEUTRALITY 

In order to determine the behavior of the anode fall, boundary conditions have to be 
applied to equations (7-9) or to set (10). In the case of the positive column and the cathode 
fall region, it is known [Refs. 45, 46, 47, 48] that for the isothermal case, the spatial 
derivatives of the particle density and the drift velocity become singular if the drift velocity 
attains the ambipolar sound speed [Ref. 49, pp. 77-86; 50], This apparent singularity appears 
near the wall, marking the end of the quasi-neutral region and the beginning of the space- 
charge sheath adjacent to the wall; the ambipolar solution breaks down because the space 
charge density cannot become larger than the local density of ions or electrons. [Ref. 50; 33; 
52: pp. 202-204], Hence, researchers have always chosen the ion sound speed as the 
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appropriate bound for the ion velocity. This condition seems a natural compromise between 
Bohm’s criterion for the collisionless sheath which requires the ion velocity exceeds or is at 
least equal to the ion sound speed and Persson’s criterion for an ambipolar flow at the plasma 
boundary, which states that the ion velocity is in fact less than the ion sound speed. However, 
these criteria are incompatible as was shown by Godyak and Sternberg. [Ref. 34] For 
example, in the numerical treatment of the plasma-wall problem without plasma-sheath 
separation, the ion velocity reaches ion sound speed in the region where plasma neutrality is 
violated [Ref. 51: p. 51]. 

It seems that since there is no strict boundary between the plasma and the sheath, it 
is not obvious which boundary conditions work the best; much controversy remains for the 
negative electrode. Attempts have been made to bridge or to remove these apparent 
discontinuities by introducing a transition layer [Ref. 27] or by choosing various plasma- 
sheath boundary conditions [Refs. 45, 34]. For example, Godyak and Sternberg [Ref. 34] 
show that the usual choice of choosing the ion sound speed as the boundary, leads to 
discontinuities of the potential, plasma density and velocity gradients. These workers 
therefore introduced as their boundary condition a specific value of the electric field that 
provides a breakdown of the neutrality at the plasma boundary and so treated the sheath and 
ambipolar model separately, attaining a smooth transition between sheath and plasma - for the 
cathode region. In their model, however, only constant ionization was assumed for species 
production. 

Fundamentally, the issue of transition in the cathode region remains somewhat 
controversial, a subject even more unclear when the positive electrode is considered, where 
essentially scant information is available. But a transition to ambipolar diffusion must occur 
[Ref 52: pp. 64-68], a transition where the electric field domination shifts to that of a density 
gradient driven region, an area where now the length scale for the number density becomes 
important. 

To capture this effect and to retain the role of charge production in the dynamical 
behavior of the system, equations (7) are again differentiated allowing the species densities 
to be written in terms of the production terms. Applying the quasi-neutral condition. 
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n~n e , then gives the transition model, equation (11) which can also be obtained from the 
sheath formulation, set (10). 
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Equation (11) 



In comparison to equations (10a- 10b), the space-charged model, the effects of 
diffusion are clearly seen not only on the rate of change of the number density or density 
gradient term (first term of (11)), but also on the energy balance terms, recombination 
(second tarn of (1 1)) and ionization (third term of 1 1). Furthermore, coulombic interactions 
are no longer evident in the energy balance, playing only a secondary role in the density 
gradient description (first term of (11)). 

A key character of equation (1 1) is that ambipolar diffusion is not yet fully formed, 
a consequence of the transition region in which perhaps the instability of the sheath, 
electrostatic repulsion of ions, is carried forth by gradient effects so that ions are further 
pushed from the positive plate. The growth in densities away from the plate then affects then- 
free diffusion, a process evolving into ambipolar diffusion. 

C. AMBIPOLAR REGION 

When the density of electrons and ions becomes large enough, their mutual Coulomb 
fields affect their free diffusion. This modification can be written in terms of the space charge 
field as given by: 
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Equation (6) 



From (6), it is seen that should the species flux densities not be equal, a significant charge 
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imbalance would arise. For example, if the plasma is much larger than a Debye length, it must 
be quasi-neutral and the diffusion rates of the species would adjust themselves so that the 
same rate from the plasma into the region of quasi-neutrality is achieved: the electrons have 
a larger thermal velocity as compared to the ions and tend to leave the plasma first. A 
positive charge is thus left behind creating an electric field of such polarity as to retard the 
loss of electrons and to accelerate the loss of ions: the plasma, by its very definition, is a 
neutral fluid. [Ref. 38: p. 139] From this formulation, the ambipolar diffusion coefficient can 



be determined as 



D = 



D '** D P' 

F, + he 



where the Einstein relation may be applied to write 



mobilities in terms of diffusions. It should be pointed out that the ratio D/p represents a 
measure of energy [Ref. 52: pp. 60-63], and that if T e =T ; , as is assumed in this work, the 
ambipolar electric field tends to enhance the diffusion of ions by a factor of two, a diffusion 
that is primarily controlled by the slower species: D a «2Dj. The primary function of the 
ambipolar field is then to serve as an energy-exchange mechanism for transferring random 
kinetic energy from the electrons to the ions such that the ion diffusion continues away from 
the positive plate. [Ref. 45] 

The model for the ambipolar region is thus formulated, again following naturally from 
equation (10) or from equations (7-9), the sheath model: proceeding as before, equations (7) 
are differentiated with respect to the spatial variable allowing set (9) to contribute to the 
densities. However, unlike previous, now set (7) is multiplied by the respective species 
diffusion coefficient, the result which gives one second order, ordinary differential equation 
for the density in terms of not only the necessary chemistry, but also in terms of the ambipolar 
diffusion, D a : 
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Equation (12) 



In comparison with equation (11), it is immediately seen that the density gradient 
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contribution, the rate of change of number density, is no longer evident. The only terms 
remaining in (12) are those of charge production; the system chemistry. As the flow 
approaches equilibrium near the plasma boundary, the number density must approach that of 
Saha so that the definition of “plasma” is upheld. 

The physical model in summary: the electric field pushes the flow from the electrode 
surface toward the density gradient dominated region in which the density gradient pushes the 
system toward the plasma boundary where eventually energy conservation must be satisfied. 
Ionization must balance recombination such that Saha equilibrium is attained at the plasma 
boundary. It is in this region that the largest of the length scales has influence, the current 
or ionization scale. This last scale defines the extent of the anode fall region, a region where 
flow chemistry becomes paramount, determining plasma equilibrium. 

But is the anode fall stable? Previous attempts [Ref. 53] at numerical analysis have 
proven challenging at best. So a better understanding of the nature of equations (7-9) and 
sets (10-12) is needed: in the next chapter, non-linear analysis techniques are applied to 
mathematical models presented previously. 
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V. ANALYSIS 



A. LOCAL ANALYSIS 

To begin the analysis of the anode fall, non-linear methods as briefly outlined in 
Appendix C are used: further information on particular ideas can be found in references 54 
and 55, among others. 

1. Space-Charge Region: Sheath 

To begin an analysis of either equations (7-9) or equations (10), non-linear methods 
as outlined in Appendix C are used: first, equations (10) are written as a system of first order 
differential equations, then stagnation points within both systems, equations (7-9) and the first 
order formulation of equations (10), are determined. Next, the local behavior about these 
stagnation regions in terms of system characteristics is evaluated, leading to a qualitative 
stability picture of the system dynamics. 

a. Equilibrium (Stagnation) Regions 

First equations (7-9) and equations (10), the anode sheath models are 
analyzed (non-dimensional parameters are given in Table 2): 
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The original model, equations (7-9) can immediately analyzed, whereas equations (10) 
have to be recast as a set of coupled, first-order differential equations. 
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Equation (10a) 

Equation (10b) 
Equation (10c) 



To find the stagnation points within these first-order system, the derivative terms are each set 
to zero and non-dimensional parameters from Table 2 are applied. The results give for the 
equilibrium criteria: 

1. From equations (8) and (10c), n j * =n e * =n* as expected; this result is then 
incorporated into equations (9a-9b) as well as (10a- 10b). 



_ _ v v 

2. Similarly, from equations (9a-9b) and (1 0a- 10b), n* = 0 or n " - — - — ; a 

a 2 n o 

stagnation point appears when n* =n* -n* = 0 and when the plasma boundary 

is approached. The latter condition immediately emphasizes the flow chemistry, 
ionization and recombination. Without recombination, the second stagnation would 
not be present, as seen from either equations (9a-9b) or equations (10a- 10b) with 
d/dy[ ]=0. The former equilibrium, on the other hand, seems artificial at first, unless 

the species densities behave such that they cross when nj =n’ =n" =0 , 

immediately implying a quasi-neutral condition since equation (8) or equation (10c) 

require n f * =n e * =n* at stagnation. In general, the fixed points found stem purely 

from the chemistry behavior, equations (9a-9b) and consequently (10a- 10b) in which 
ionization and two-body recombination are the only processes included. If, on the 
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other hand, three-body recombination is added to the equations, then additional 
stagnation regions become evident since equations (9a-9b) are now cubic. The 

equilibrium region near n’ =« ’ = 0 still remains, however. Similarly, if 

ionization is the only production term considered in the sheath model, then the 
stagnation region near the plasma boundary “disappears.” The condition 

n g * =n. m =n* = 0 remains for all cases: ionization only, ionization with two-body 

recombination, ionization with two- and three-body recombination. So it appears that 
continuity (equations (9)) requires flow equilibrium when perhaps quasi-neutral 
conditions are initially met, at the sheath boundary, and when the plasma interface is 
approached. It is recombination that establishes the latter while the presence of 
ionization stipulates the former. 



3. The current densities at the equilibrium position can be found by using n ’ = 0 



v v _ v v 

or w* = with equations (7a- 7b): if «’ = , then 
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since n -n -n when n' -n' -n' Here, the influence of not only 

o t o t o i e ? j 

ionization and recombination, but also of diflusion on the current density distributions 
are clearly seen. Similarly, from equations (7a-7b), j i e * =0 V« ’ =0 



4. Charge conservation is used to isolate the value of the electric field at stagnation 

v o v 

near the plasma boundary, where n* = . That is, 

a 2 n o 



V-J=0 =>J=j e+ j=J o or j e +J,= \ gives 
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from which E * is obtained graphically since the right-hand side of this expression 
simply represents a constant. Applying the Arrhenius function [Ref. 52: pp. 94-109] 



for the ionization. 




and using numerical data for nitrogen 



yields E*~ 0.54 . From this, y' e *®0.99816, j. *=<0.0005445, w* = 0. 15695 

These results show the relative magnitude of the electron and ion currents within the 
anode fall, j e »j ; as expected. If, on the other hand, j . e * = 0 V n ’ = 0 is used, an 

inconclusive result is obtained: E * cannot be isolated from either equations (7-9) 



or equations (10). Further, 




does not permit zero field. Therefore, 



it seems reasonable that a finite field should exist over the entire range of the anode 
sheath. Von Engel [Ref 56: pp. 12-20] highlights such arguments: a finite field must 
be allowed, which could then introduce a potential trough. A virtual anode is thought 
to exist if the space charges have the same sign as the nearby electrode, with the 
emitting surface located at the potential minimum [Ref. 56: pp. 14-17], 



The fixed points found in the model represented by equations (7-9) and in part by set (10) are 
listed in Table 5 below, featuring the role of ionization, recombination and diffusion in the 
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existence of the fixed point near the plasma boundary. Again, continuity (equations (9)) 
seems to require flow equilibrium when perhaps quasi-neutral conditions are initially met, the 
sheath boundary, and when the plasma interface is reached. It is the recombination 
phenomena that establishes the latter stagnation region while the presence of ionization 
stipulates the former. 
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Table 5: Stagnation Regions in General Sheath Model - Ionization with Two-body 

Recombination 
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b. Jacobian at the Equilibrium Positions 

To determine the system characteristics about its fixed points, a local Taylor 
series expansion is carried out about the equilibrium of interest (Appendix C); the rate of 
change of the flow is given in terms of the Jacobian, which for equations (7-9) becomes: 
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Here, six variables affect the system dynamics: ion and electron number densities, 
electric field, species diffusivities, and ionization. Of these, the diffusion coefficients are 
assumed constants, so essentially four variables drive the behavior: species densities, electric 

field and consequently ionization. If A jnner is evaluated at the stagnation points, the 

resulting eigenvalues describe the dynamics in the neighborhood of the fixed point (Appendix 
C). Specifically, at the plasma boundary equilibrium with nitrogen data: 

A, = 1.376704 A 2 = 0 

X 3 = 0.630455 X 4 = -1.37437 

X 5 = - 0.63279 

For this case, the eigenvalues are real and distinct, leading to a multi-dimensional type of 
“saddle” behavior (Appendix C). The stagnation region near the plasma boundary is unstable: 
the largest eigenvalue is positive, indicative of exponential growth. Appendix D lists the roots 
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of the characteristic polynomial of A jnner . It appears A, gives exponential growth - the 



process of ion repulsion from the positive plate forces the instability of the sheath. The sheath 
drives the system toward the first stagnation criterion, toward quasi-neutrality. 

Provided a finite field exists when n e * = n* = n * = 0 , then the fixed point in its 



vicinity gives a similar qualitative behavior: the eigenvalues of the resulting Jacobian are given 
by 
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These eigenvalues are always real with two characteristics equal while the remaining are 
distinct: Table 6 lists some numerical examples for nitrogen data. 
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Table 6: Eigenvalues of Stagnation Region n,=n e =n=0 with Finite Field - Numerical Values 

for Nitrogen 

Again, the essential feature of the sheath is instability: ions are driven from the 
positive wall as required by electrostatics. The stagnation region near n*=n*=n*= 0 is 

unstable, perhaps pushing the system away toward the recombination induced equilibrium 
near the plasma boundary. It is the strength of the initially applied electric field that imposes 

the stability characteristics for the fixed point n e * =n*=n* = 0 . This leads to the 

following question: is the sheath model, as presented by either equations (7-9) or set (10), 
an example of a system in which a certain type of “fixed point” is always present, a stagnation 
region satisfying quasi-neutrality and thereby representing the sheath boundary? If so, does 
this equilibrium present a type of transcritical bifurcation (see Appendix C) whose stability 
perhaps depends on the initially applied electric field? The foregoing does suggest so. 

The precarious nature of the fixed points found in the sheath model suggest that 
numerical integration methods will not prove fruitful over the entire domain of the anode fall. 
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However, by exploring the various regimes of equations (7-9) in the form of the second order 
differential equations given by equations (10), further insight into the anode fall could be 
expected. 

2, Transition Zone: Quasi-Neutrality 

Previously, it has been shown that the sheath is inherently unstable, that the ions are 
electrostatically repelled from the electrode (anode) surface. If through that process quasi- 
neutrality occurs, then for steady-state conditions either density gradients or diffusive 
processes have to drive the system toward the plasma boundary since the condition of 

n~n e stipulates the electric field be constant; so something else has to continue driving 

the ions away from the anode wall. Equation (11) highlights the density gradient and is 
considered next as a set of coupled first-order differential equations: 
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Equation (11) 



where the electric field is constant due to quasi-neutral conditions, 



33 



a. Equilibrium (Stagnation) Regions 

To find the stagnation points within these first-order system, the derivative 
terms are each set to zero and non-dimensional parameters from Table 2 are applied. The 
results gives for the equilibrium criteria: 

1. Setting d [ ]/dy=0 in equation (11) immediately gives the previously found 

v v 

stagnation condition near the plasma boundary: n * = The fixed 

a 2 n c 

point n * = 0 again seems to present a mathematical phenomenon since equation 

(11) does not describe the sheath condition. If quasi-neutral arguments are again 

used, then perhaps the condition n * = 0 may have meaning - the plasma shields 

itself against external influences (electric field), and so, quasi-neutral characteristics 

are fundamental to its definition. The stagnation region at n ’ = 0 does, perhaps, 

present a physical phenomenon: the onset of quasi-neutrality or the sheath boundary, 
found as a consequence of ionization present in the sheath model and in equations 
(11), the “transition” model. 

2. The magnitude of the constant electric field, and therefore constant ionization, at 
the start of the quasi-neutral transition is determined by the matching conditions when 

n~n e . These criteria are either found through analytically matching equations ( 1 0) 

with equations (1 1) for some small parameter(s), or numerically. Numeric methods 
will be considered later in this work. Suggestions toward an analytic matching 
attempt are provided in paragraph B below. 

b. Jacobian at the Equilibrium Positions 

To determine the system characteristics about its fixed points, a local Taylor 
series expansion is carried out about the equilibrium of interest (Appendix C); the rate of 
change of the flow is given in terms of the Jacobian, which for equation (11) becomes: 
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Since the electric field becomes constant at the onset of quasi-neutrality at the sheath 
boundary, the electric field thus locks in the behavior of the transition region. To determine 

these stability characteristics, the general eigenvalues of A tra . . are determined: 



x _ aE(D-D t ) ^ ) ja 2 E 2 (D r D e ) 2 -4v(D i+ D e Kbi+dg) + %HD i +D t )(bj * dh) 

U " 2{D i+ D e ) 2(D t+ D e ) 

where the non-dimensional parameters are listed in Tables 2 and 10. These eigenvalues 
determine the dynamical behavior of the transition when A transition is evaluated at a 

fixed points, albeit the electric field strength is determined from the sheath as the system 
reaches n~n e . By inspection of the first term found in these eigenvalues, the following 

is observed: 

1 . If the field is positive, then the system displays exponential decay, stability; but 
electrostatics demands that ions be repelled from the positive plate. This leads to the 
expectation that the transition zone is also unstable, similar to the sheath. 
Consequently, the field should be negative at the point of quasi-neutrality. 

2. Using the stagnation criteria for n ' = 0 with finite field indicates that the 
characteristics should display stable behavior for positive field, the diffusion difference 
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(D r D e ) is negative, and unstable behavior for negative field. The initially applied 
electric field strength then does determine the stability of the fixed point n * = 0 , 
an equilibrium which appears when ionization is present in the model. 

- V 

3. To find the behavior near the plasma boundary when n ’ = , nitrogen 

a, « 

2 O 

data along with an arbitrary selected field strength, £(_ y-_- ) = -5.00 are used to 

give the eigenvalues Aj—0.041 and A 2 =5.026: the eigenvalues are real, distinct and 
predict a saddle-type behavior near the plasma boundary - a robust linearization (see 
Appendix C). For this case, the flow field, locally, is characterized by a very strong 
unstable manifold corresponding to A 2 . A similar result was found by analyzing the 
sheath model presented by either equations (7-9) or equations (10). 

4. On the other hand, if the sheath behavior at quasi-neutral formation gives a 
positive field, e.g., E(y- ) = 5.00 , then the Jacobian eigenvalues predict stability 

locally about the fixed point: ^=0.021 and A2=-5.013 where now the saddle behavior 
is marked with a strong, stable manifold, contrary to the sheath analysis earlier. 

Thus, the criteria found at the onset of quasi-neutrality at the sheath boundary form 
the initial conditions imposed upon the transition zone and hence govern the stability 

- _ v v 

characteristics of the fixed points at n * = 0 and n * = — - — 

a,» 

2 o 



3. Ambipolar Region 

In the sheath, ions are repelled electrostatically until quasi-neutrality occurs, where 
chemical processes balance such that the ion density approaches that of the electrons - the 
electric field becomes constant. If this balance does not reflect total charge neutrality, i.e., 
Saha equilibrium, then other processes (the field is now constant) must drive the system 
toward such equilibrium; however, a stabilizing mechanism must be available once the 
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